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focusing  model  for  the  ions  is  presented. 
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TIME  DEPENDENT  ELECTRON  AND  ION 
FLOW  IN  PINCHED-BEAM  DIODES* 


INTRODUCTION 

In  order  to  gain  a better  understanding  of  the  dynamic  processes  taking  place  in  high- 
aspect-ratio  pinched-beam  diodes,  a two-dimensional  computer  simulation  code  DIODE2D  1 has 
been  developed  to  study  these  processes  for  a wide  range  of  diode  parameters.  This  code, 
which  follows  the  relativistic  dynamics  of  electrons  and  ions  coupled  to  their  self-consistent 
electric  and  magnetic  fields,  contains  a fully  time-dependent  emission  scheme  which  allows 
time  dependent  diode  phenomena  to  be  studied.  Figure  1 shows  the  diode  geometry  used  in 
the  simulations.  Both  flat  faced  and  hollow  cathodes  were  simulated.  For  flat  faced  cathodes 
(Fig.  la)  electrons  are  emitted  self-consistently  all  along  the  whole  cathode  (K)  while  for  hol- 
low cathodes  (Fig.  lb)  electrons  are  emitted  only  over  a region  of  radial  width  SR.  Since  the 
version  of  the  code  used  in  these  studies  allowed  only  a rectangular  simulation  region,  the  po- 
tential along  the  remaining  segments  of  the  boundary  was  graded  to  emulate  the  actual 
geometry.  Emission  of  the  ions  from  the  anode  (A),  also  self-consistent,  was  usually  over  its 
full  face  out  to  the  outer  radius  of  the  cathode. 

In  this  report,  we  first  look  at  a composite  picture  of  the  impedance  of  the  diode  which  is 
the  result  of  many  simulations.  This  description  of  the  impedance  is  then  used  to  analyze  ex- 
perimental data  from  two  representative  shots,  one  taken  on  Gamble  I and  the  other  taken  on 
Gamble  II.  The  diode  current  for  the  shots  is  calculated  as  a function  of  time  using  the  experi- 
mentally measured  corrected  voltage  as  an  input  parameter.  The  calculated  currents  are  then 
Manuscript  submitted  January  18,  1978. 
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compared  to  the  measured  currents  for  the  two  shots. 

The  dynamic  details  of  one  of  the  simulations  used  to  construct  the  impedance  picture  is 
then  presented.  The  collapsing  hollow  ring  and  pinch  formation  are  shown,  as  well  as,  the  to- 
tal diode  current  and  ion  current  as  a function  of  time.  We  then  model  and  show  the  project- 
ed ion  focus. 

DIODE  IMPEDANCE: 

In  our  numerical  studies  the  primary  goal  at  this  stage  was  to  determine  the  impedance  of 
the  diode  as  a function  of  various  diode  parameters.  Limitations  of  computational  speed  make 
it  impractical  to  stimulate  the  complete  time  history  (50-100  nsec)  of  the  pinched-beam  diode. 
Simulations  describing  5 nsec  of  real  time  can  be  calculated,  however,  and  this  is  sufficient 
time  for  an  ion  induced  pinch  to  form  for  diodes  with  radii  less  than  approximately  10  cm  if 
the  dynamics  of  cathode  and  anode  plasma  generation  are  ignored.  We  will  assume  that  at  the 
beginning  of  the  simulation  both  the  cathode  and  anode  plasma  exist  over  all  or  some  part  of 
the  cathode  or  anode.  This  control  of  the  actual  emission  regions  is  used  to  model  the  regions 
where  actual  plasma  generation  and  electron/ion  emission  take  place.  In  the  emission  regions 
the  electrons  and  ions  are  emitted  in  a time  dependent,  space-charge-limited  manner.  The 
electron  and  ion  dynamics  are  then  followed  until  a quasi-steady  state  is  reached,  if  ever  It 
should  be  noted  that  in  some  of  the  data  presented  here  the  simulations  did  not  reach  a 
steady-state  although  the  total  current  reached  a quasi-steady  value. 

Two  types  of  diodes  were  studied.  The  first  has  a full,  flat  faced  cathode  with  a maximum 
radius,  R,  of  5 cm  and  R/D's  from  4 to  30,  where  D is  the  gap  spacing.  The  second  is  a very 
hollow,  flat  faced  cathode  with  R of  6.4  cm,  SR/R  of  0.1 1,  and  R/D  of  20.  Here  A R is  the  radi- 
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al  width  of  the  emitting  surface  on  the  cathode.  Both  diodes  were  simulated  for  cases  with 
both  electron  and  ion  currents  and  then  electron  current  only.  The  total  current  for  one  diode 
in  which  ions  were  emitted  only  from  R =0  to  R = 3/4 R was  also  calculated  at  an  R/ D of 
10.  Figure  2 shows  a plot  of  the  total  diode  current  scaled  by  D/R  as  a function  of  voltage  for 
each  of  these  diodes.  The  dashed  lines  connect  results  obtained  at  an  R/D  of  20.  The  closed 
and  open  squares  represent  data  for  a solid  cathode  with  and  without  ion  current.  The  shaded 
region  connecting  triangular  points  represents  the  computed  diode  current  for  a hollow  cathode 
with  1R/R  o/O.ll.  The  upper  bound  in  this  case  being  for  electron  plus  ion  current  while  the 
lower  bound  in  this  case  being  for  electron  plus  ion  current  while  the  lower  bound  is  that  for 
electron  current  only.  The  plus  symbol  represents  the  calculation  for  the  case  where  the  ion 
emission  was  only  for  radii  smaller  than  3/4  R.  Also  plotted  for  reference  in  Fig.  2 are  the  sa- 
turated parapotential  current  8.5  (kA)  y In  [y  + (y2  — 1)  1/2 ] RID  (curve  PP)  and  the  criti- 
cal current  for  self  pinch  8.5  (kA ) (y  2 — 1 ) 1//2  R/D  (curve  CC)  curves. 

In  the  upper  half  of  Fig.  3 the  experimental  voltage  and  current  curves  are  shown  for 
shot  5707  on  Gamble  I.  The  very  hollow  cathode  used  had  an  outer  radius  of  42  mm  and  a ra- 
dial emission  width  of  7 mm.  After  whisker  exposion  the  corrected  gap  was  taken  to  be  3 mm 
and  an  average  plasma  closure  velocity  of  1.8  cm//usec  was  inferred  from  an  analysis  using  the 
so  called  ''semi-emperical”  current  formula. 2 The  lower  half  of  Fig.  3 shows  the  effective  gap 
and  the  measured  impedance  as  a function  of  time.  The  initial  plasma  motion  is  assumed  to  be 
correlated  with  the  onset  of  the  current  and  hence  delayed  somewhat  from  the  voltage  rise. 
The  gap,  as  a function  of  time,  is  given  by  D (T)  =3.0  — 0.018  (T  — 9 ),  where  D is  in  mil- 
limeters and  7"  is  in  nanoseconds. 

Figure  4 shows  the  voltage,  current,  and  impedance  curves  for  shot  1776  on  Gamble  II. 
Another  very  hollow  cathode  was  used  with  an  outer  radius  of  65  mm  and  a radial  emission 


LEE  AND  GOLDSTEIN 


width  of  4 mm.  The  initial  corrected  gap  was  taken  to  be  4.4  mm  and  the  average  plasma  clo- 
sure velocity,  calculated  from  the  time  to  short  circuit,  was  4.0  cm//asec.  The  gap,  as  a function 
of  time,  is  then  given  by  D ( D =4.4  — 0 .04  (7"  — 10). 

Table  l gives  a summary  of  the  important  parameters  for  the  two  shots.  Shot  5707  was 
analyzed  for  times  from  40  nsec  to  100  nsec  while  shot  1776  was  analyzed  from  40  nsec  to  100 
nsec.  77,  is  the  time  at  which  the  pinch  was  formed  on  axis.  The  effective  R/D  varied  with 
time  as  the  gap  closed.  In  shot  5707  R/D  was  20  at  40  nsec  and  became  infinite  when  the 
diode  sorted  out.  Since  all  the  simulations  for  hollow  cathodes  were  at  R/D  = 20  it  is  implicit- 
ly assumed  that  to  lowest  order  the  current,  at  a fixed  voltage,  increases  linearly  with  R/D. 
This  seems  to  be  true  for  the  very  hollow  cathodes. 

In  Fig.  5 the  normalized  current  versus  the  corrected  voltage  is  plotted  for  the  two  shots. 
The  upper  dashed  line  is  the  total  current,  electron  plus  ion,  for  the  R/D  = 20  very  hollow 
cathode  simulations  given  in  Fig.  2.  The  lower  dashed  line  is  the  electron  only  current  from 
the  same  simulations.  The  two  curves  plotted  on  the  upper  part  of  the  figure  are  the  evolution 
in  time  of  the  normalized  current  versus  voltage  for  the  Gamble  I and  Gamble  II  shots.  The 
arrows  along  the  curves  point  in  the  direction  of  the  advancement  of  time.  The  vertical  arrows 
indicate  the  time  at  which  the  pinch  arrived  on  axis.  In  the  analysis  of  the  data  the  linear  gap 
closure  model  discussed  previously  was  used. 

Since  ion  flow  in  the  diode  is  necessary  for  pinch  formation/  one  would  expect  full  elec- 
tron and  ion  flow  to  have  developed  at  about  the  time  the  pinch  forms  on  axis.  The  voltage- 
current  characteristics  would  then  be  expected  to  follow  the  upper  dashed  curve  on  the  figure. 
Shot  5707,  the  Gamble  I shot,  follows  the  simulation  curve  remarkable  well  shortly  after  the 
pinch  forms  on  axis.  The  data  for  the  Gamble  II  shot,  1776,  although  in  the  general  vacinity  is 
somewhat  lower  than  the  predicted  value. 
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The  key  feature  in  the  analysis  of  the  data  is  the  model  for  the  gap  closure.  The  experi- 
mental data  cannot  be  analyzed  directly  since  the  gap,  as  a function  of  time,  cannot  be  meas- 
ured directly.  In  order  to  compare  the  experimental  data  with  the  simulations  we  must 
somehow  model  the  dynamics  of  the  gap.  Of  all  the  areas  of  diode  physics  the  formation  of 

the  cathode  plasma  by  whisker  explosion  and  its  subsequent  motion  is  probably  the  least  un- 
derstood. It  is  usually  assumed  that  in  the  first  few  nanoseconds  of  the  pulse  the  field  emis- 
sion current  flowing  in  the  cathode  whiskers  heats  them  until  they  explode  forming  a cathode 
plasma  sheath.  The  formation  of  the  sheath  and  its  resultant  expansion  velocity  seems  to 
depend  critically  on  the  initial  rise-time  of  the  electric  field. 2,4  For  lack  of  any  better  idea  one 
usually  trys  to  find  the  average  expansion  (closure)  velocity.  This  seems  to  work  fairly  well 
with  the  lower  power  machines  such  as  Gamble  I and,  hence,  the  good  agreement  in  Fig.  5. 

On  Gamble  II,  however,  there  are  other  competing  mechanisms:  namely,  the  heating  of 
the  cathode  by  electrons  returned  by  the  time  dependent  electromagnetic  fields  due  to  induc- 
tive effects  in  the  diode,  and  from  ion  current  striking  the  cathode  with  current  density  in- 
versely related  to  the  diode  gap  and  radius.  Both  of  the  mechanisms  could  increase  the  closure 
velocity  as  a function  of  time.  In  light  of  these  effects,  we  have  investigated  a simple  model 
where  we  arbitrary  assume  that  the  plasma  velocity  increases  linearly  in  time  from  some  initial 
value.  The  difference  between  the  average  value,  which  can  be  calculated  from  the  time  of  the 
short  circuit,  and  the  initial  value  is  a parameter  in  this  model.  Figure  6 shows  a comparison 
between  the  constant  velocity  model  and  the  linearly  increasing  velocity  model.  The  difference 
between  the  two  models,  since  both  have  the  same  average  velocity,  is  that  the  gap  closes 
more  slowly  at  first  in  the  latter.  Using  a value  of  a equal  to  0.5,  we  have  re-analyzed  shot 
1776  and  the  results  are  seen  in  the  lower  figure  on  Fig.  5.  Using  this  model  then  gives  a 
much  better  fit  of  the  experimental  data  to  the  simulation  results.  This,  of  course,  in  no  way 
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proves  that  the  gap  closes  in  a nonlinear  fashion  but  since  there  is  other  evidence  for  the  pre- 
viously discussed  mechanisms  this  just  re-enforces  ones  belief  that  plasma  gap  closure  is  a 
complicated  and  not  very  well  understood  phenomenon. 

TIME  DEPENDENT  PINCH  FORMATION  AND 
ION  GENERATION: 

In  this  section  details  of  the  time  dependent  formation  of  a pinch  from  a very  hollow 
cathode  are  presented.  The  parameters  for  the  simulation  are  as  follows:  the  diode  voltage, 
750  kV,  the  outer  radius  of  the  cathode,  6.4  cm,  the  radial  width  of  the  emitting  surface,  7 mm, 
and  the  anode-cathode  gap  was  3.2  mm.  Ions  were  emitted  over  the  whole  anode  out  to  the 
maximum  cathode  radius.  Figure  7 shows  the  diode  geometery  and  a composite  picture  of  the 
electron  current  density  collected  on  the  anode  plane  as  a function  of  time.  The  collapsing 

hollow  ring,  which  reaches  the  axis  and  forms  the  pinch  in  about  2.5  nanoseconds,  is  clearly 
seen.  Although  the  peak  current  density  in  the  pinch  region  is  extremely  high,  of  the  order  of 

1.6  MA/cm  2,  only  about  40  kA  of  total  electron  current  are  contained  within  the  pinch.  The 
pinch  itself  was  not  stable  in  this  simulation  but  would  disappear  soon  after  forming.  Then  at 
some  intermediate  radii,  between  2.5  - 4.5  cm,  one  or  two  hollow  rings  would  reappear.  They 
would  then  collapse  in  about  1.5  to  2 nanoseconds  to  re-form  the  pinch. 

The  small  peak  in  current  density  at  large  radii  is  most  probably  due  to  electrons  emitted 
from  the  inside  edge  of  the  cathode.  The  maximum  of  this  peak  is  a current  of  about  15 
kAlcm  2 and  amounts  to  a current  of  about  60  kA,  or  more  than  that  in  the  pinch  itself.  The 
remaining  current,  approximately  250  - 300  kA  is  spread  out  more  or  less  uniformly  in  radius 
over  the  anode. 


6 


f NRL  MEMORANDUM  REPORT  . ‘02 

The  total  diode  current  and  the  ion  current  collected  on  the  cathode  are  shown  in  Fig.  8 
as  a function  of  time.  The  first  thing  to  take  note  of  is  that  the  total  diode  current  remained 
very  steady  for  the  duration  of  the  simulation  at  a value  of  about  550  kA.  On  the  other  hand 
the  ion  current,  and  hence  the  electron  current  also,  oscillated  in  time.  The  average  value  of 
the  ion  current  was  around  140  kA  giving  an  average  ion  current  to  electron  current  ratio  of 
1/3.  The  ion  current,  which  would  peak  some  0.2  - 0.3  nanoseconds  before  the  pinch  formed 
on  axis,  had  maximum  values  near  170  kA. 

In  Fig.  9 the  density  of  ion  current  collected  by  the  cathode,  at  the  time  of  the  first  max- 
imum of  the  total  ion  current,  is  plotted  as  a function  of  radius.  This  figure  is  representative  of 
the  ion  current  density  during  the  collapse  of  the  e'ectron  hollow  ring  on  the  anode.  The  peak 
of  the  ion  current  density,  which  has  a value  of  about  9 kA/cm2,  is  located  at  the  radial  posi- 
tion that  the  electron  ring  had  on  the  anode  when  the  ions  were  emitted,  some  0.4 
nanoseconds  earlier.  The  average  ion  current  density  over  the  whole  cathode  is  — 2 kA/cm  -. 

We  now  use  the  results  of  this  simulation  to  model  the  ion  focusing  by  a shallow  spheri- 
cal anode  with  a radius  of  curvature  of  12.5  cm.  The  planar  geometry  used  in  the  simulation 
was  mapped  into  the  spherical  dish  geometry  used  in  the  acual  experiments. 5 The  ion  orbits 
calculated  by  the  code  were  then  geometrically  projected  into  the  region  outside  of  the  non- 
emitting portion  or  the  cathode.  The  variation  of  ion  directions  at  the  cathode,  about  2 de- 
grees, were  taken  into  account  when  the  projections  were  calculated  and  the  ion  current  densi- 
ty was  averaged  over  a complete  oscillation  period.  Figure  10  now  shows  the  axial  variation  of 
the  radial  current  density  profiles  as  one  moves  to  either  side  of  the  anticipated  focal  position. 
At  the  focus  the  peak  current  density  of  ^60Avl/cm  2 is  over  a radius  of  about  4 mm.  Also, 
notice  that  while  the  current  density  at  large  radii  is  small  it  contains  about  75%  of  the  total  ion 
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current  (120  kA  time  averaged).  Some  of  this  current  could  be  brought  into  focus  by  using  as- 
pheric surfaces.  However,  time  dependent  fluctuations  will  probably  put  a limit  on  the  smallest 
focal  size  attainable. 

CONCLUSION: 

We  have  shown  that  the  time  dependent  impedance  of  the  diode  can  be  reasonably  inter- 
preted using  the  results  of  the  computer  simulations.  And,  that  pinch  formation  and  other 
dynamic  processes  that  take  place  in  relativistic  pinch-beam  diodes  can  be  studied  in  detail. 
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Fig.  1 — Diode  geometry  for  simulations  Fig.  2 — Total  diode  current  normalized  by  D/R 

as  a function  of  diode  voltage.  The  parapoiential 
current  (PP>  and  critical  current  (CC)  are  also 
shown  for  comparison 
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Fig.  3 — Corrected  diode  voltage  and  current  for  shot 
5707  In  lower  graph  the  experimental  value  of  the 
impedance  is  shown  along  with  the  gap  closure  used 
to  reduce  the  data.  • 
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Fig.  4 — Corrected  diode  voltage,  current,  and 
impedance  versus  time  for  shot  1776. 
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Fig.  5 — Normalized  current  calculated 
from  the  computer  simulations  versus  vol- 
tage for  the  two  shots.  The  upper  graph 
uses  the  linear  gap-  closure  model  and  the 
lower  graph  uses  the  nonlinear  model. 


11 


NRL  MEMORANDUM  REPORT  3702 


LEE  AND  GOLDSTEIN 


i 

\ 


30  r 
60}- 
40)- 


Fig.  10  — Radial  proton  current  density  profiles  from 
time  averaged  computer  simulation. 
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